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■ Abstract. We compare different driving scenarios for the control of 
' matter-wave squeezing of a coherent preparation in the vicinity of the 

■ hyperbolic instability of the two-mode bose-Hubbard model. In par- 
C/3 ■ ticular we contrast the quantum Zeno suppression of decoherence by 

broad-band erratic or noisy driving, with the Kapitza effect obtained 
for high frequency periodic monochromatic driving. 

» ■ 

^: 

^3 '■ 1 Introduction 

Si 

. I The physics of Bosc-Einstcin condensates (BECs) confined in periodic lattice poten- 

. tials is captured in the tight binding approximation, by the Bose-Hubbard Hamilto- 

enian (BHH) |1I2I3I4I5|5] . With only two coupled condensates, this model reduces to 

I , the Bose-Josephson form, aka the Bose-Hubbard dimcr, written in spin form as 

pH ■ H = -[K + fit)] J, + C/,/f , (1) 

, where K is the inter-mode coupling, which we later modulate by adding time depen- 

'— dent driving /(t), and U is the intra- mode interaction parameter. The spin operators 
I , are defined in terms of the bosonic creation and annihilation operators a| and Oi of 

^ ■ particles in mode i = 1,2. Namely, Jz = {fii — ri2)/2 corresponds to the occupation 

2^ ■ difference, while = {a\d2 + a\ai)/2, and Jy ~ {a\a2 ~ alai) / {2i) correspond to the 

real and imaginary parts of the inter-mode coherence. Since the total particle number 
is conserved at fii + n2 = N, the total spin is set to j = N/2. For simplicity, we 
assume repulsive interaction U > 0, noting that results for an attractive interaction 
are easily deduced by taking K — —K and E — >■ —E. 

The characteristic interaction parameter of the Hamiltonian ([T]) is. 



o 

(N 



u = NU/K . (2) 



Single particle coherence, the hallmark of a global macroscopic order, is determined 
■ by the length of the Bloch vector S = (J)/j. The classical limit, which preserves 

^ , S' = 1, is obtained by restricting the quantum Hilbert space to the subset of spin 

coherent states, 

\0,(p) = cxp{-iLpX) cxp{-i0Jy) (3) 

where n = cos is the relative population imbalance, and ip the relative phase between 
the two condensates. To the extent that an initial spin coherent state evolves only 
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Fig. 1. (a) Phase-space structure of the BHH for u — 2.5. Lines depict equal energy contours. 
A separatrix with an isolated hyperbolic point at S,r = |7r/2,7r) separates the "islands" and 
"sea" regions. The Wigner distribution for the coherent state is shown in (b), and for a 
squeezed state in (c). 



to other coherent states, the dynamics can be described by the mean-field equations, 
where the spin operators in Eq. ([T]) arc replaced by c- numbers with 0(1 /TV) accuracy. 

The equal energy contours of the classical dimcr Hamiltonian in the Josephson 
interaction regime, 1 < li < are illustrated in Fig. [T] As seen, the spherical 
phase-space consists of three distinct regions of motion: Two nonlinear "islands" are 
separated from a "sea region" by a separatrix. The sea trajectories correspond to 
Rabi- Josephson oscillations of the population [11718) . whereas the island trajectories 
correspond to self-trapped motion |2l9ll0j . 

Linearization of the dimer BHH around the elliptic fixed-point So = (1,0,0), 
corresponding to the spin coherent state |7r/2, 0), gives the Josephson frequency loj = 
K\/u + 1 of small oscillations around the ground state. By contrast, linearization 
around the hyperbolic fixed-point 8,^ = (—1,0,0), corresponding to the spin coherent 
state |7r/2,7r), gives the Josephson squeezing rate 

wj = K\/u - 1 (4) 

Note that in the Josephson regime u > 1 implies the appearance of the separatrix 
and hence the hyperbolic behavior of S^. We also assume that u <C N'^, so that the 
"sea region" and the hyperbolic structure are well resolved. 

Our goal is to study the effect of driving on the squeezing of the coherent state 
|7r/2, tt), whose phase-space Wigner distribution is approximately a minimal Gaussian 
centered at S^, with uncertainty radius rg = (2/iV)^/^. Due to the hyperbolic dynam- 
ics in the absence of driving, this Gaussian undergoes expansion and contraction with 
factors e*"'-'*, respectively, along two principal axes (Fig. [T]d,c). The angle between 
the principal axes is twice the value of = tan~^ {wj/K). 

Due to the growth of variance along the expanding axis of the squeezed Gaussian 
state, the one-particle coherence S decreases. Using a simple phase space picture, a 
good approximation for the loss of coherence is given by [35] , 

S{t) = exp(^-i((r2)-r2)^ = cxp [- co^ {20) sinh^ {w jt)} . (5) 

where (r^) is the angular spreading of the squeezed Gaussian. One implication of 
this squeezing process is the hyperbolic growth of deviations from mean-field theory 
[516116] . Due to the finite phase-space volume, the short time squeezing is then fol- 
lowed by a series of revivals and collapses due to the repeated folding and interference 
of the Wigner distribution [11I12I13I14I15TT6] . 
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2 Driving scenarios 

In this work wc compare three strategies for the stabilization of 8,^ coherent prepara- 
tion by introducing a time-dependent driving f{t) as in Eq.(IT|), i.e. by the modulation 
of the two-mode coupling strength. The simplest driving scheme assume a periodic 
mono-chromatic sinusoidal time dependence: 

f{t) — Dsm{f2t + <j>), [harmonic driving] (6) 

with intensity D and frequency f2. Defining dimensionless frequency fi = i7/wj, 
and dimensionless driving strength q = ^/uD/f], fast and slow driving correspond to 
i? ^ 1 and i7 ^ 1, respectively, whereas q ^ SI and q Q correspond to strong 
and weak driving. 

Next we consider a broad-band erratic driving field /(i), with zero average and 
a short correlation time. This can be viewed as a realization of a stochastic process, 
such that upon averaging 

(fit) fit')) - 2DS{t - 1'), [erratic driving] (7) 

We stress that 'erratic driving' still refers to Hamiltonian evolution with a determinis- 
tic, well controlled f(t). An experimentalist can reproduce the same f{t) many times, 
and perform a quantum measurement of the outcome. Optionally, he can carry out 
experiments with different realizations of f{t), and accumulate statistics. 

Finally, we investigate the case of full-fledged quantum noise, induced by the 
coupling of the system to a "bath". In this case averaging is not at the courtesy 
of the experimentalist, but an essential ingredient in the proper description of the 
reduced dynamics. Assuming that the quantum noise has the same correlation as 
postulated in Eq.(l7]), the probability matrix would obey the master equation, 

P = -i[H, p] - D[X, [X, p]], [noisy driving] (8) 

where the bath degrees of freedom have been traced out. This is equivalent to aver- 
aging the Hamiltonian erratic driving dynamics over all possible realizations of f{t). 

3 Harmonic driving - Kapitza effect 

In order to understand the effect of harmonic driving, we note that for small popula- 
tion imbalance n ^ N (i.e. when 6 « 7r/2) we have Jx ~ {N/2) cosip, and the BHH 
takes the form of a pendulum Hamiltonian, 

H = Un^ ~ {1/4)N K cos ip , (9) 

The stable Sq and the unstable S^r fixed points corresponds to the 'down' and 'up' 
directions of the pendulum. Within this approximation, off-resonant high-frequency 
weak drive, corresponding to vertical driving of the pendulum axis, results in an 
effective potential = {l/A)q'^Kj siB.^{(p) }17I18I19I20] . Sufladently strong driving 
{q^ > 2) stabilizes the ip = ir fixed-point, producing the Kapitza inverted pendulum 
effect. 

This analysis was extended to the full spherical phase space and general driving 
term [5T]. Given that the Hamiltonian is of the form H + f{t)W, the effective po- 
tential is — [1/(4J7^)] [W, [W,H]], which for W = Jx contributes a term proportional 
to sin^ 9 sin (p, as well as additional terms that slightly renormalize the bare values 
of K and U . This effective potential leads to the stabilization of the Sjr point, and 
correspondingly, to the suppression of quantum squeezing, as shown in Fig. [2] 
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Fig. 2. The classical trajectories (a,c) and a representative Husimi distribution (b,d) for the 
dynamics that is generated by the BHH with high frequency harmonic driving, assuming S-n- 
coherent preparation. The parameters are u = 100, — 30, and A'^ — 100. The strength of 
the driving is g = in (a,b) and for the lower line in (e). Otherwise it is g = 3. 



4 Erratic driving - Squeezing axis randomization 



Fig.[3]shows the dynamics of the single particle coherence under the influence of erratic 
driving. In contrast to the drive- free dynamics (O, the orientation of the wavepacket 
relative to the principal squeezing axes is not constant: the erratic driving randomizes 
this orientation on time scale to = 1/{2D), resulting in angular diffusion [52]. The 
squeezing rate is no longer constant, thus Eq. ([5]) should be replaced by the general 
expression 

S{t) = exp{-rl sinh2(yl)} , (10) 



with a squeezing parameter A that accumulates in a stochastic fashion. The accumu- 
lated squeezing parameter yl in a given realization of f{t) is a sum of uncorrelated 
variables, and therefore, according to the central limit theorem, its many-realization 
distribution is normal, with some mean value fi and dispersion a. It follows from 
Eq. ([TUI that the corresponding distribution of S is log- wide, with typical value (me- 
dian) that might be very different from the mean: 



-^median 



cxp • 
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Smcan ~ CXp ^ ^ 



(P)} 

2'"' cosh(2^i) - 1 



(11) 
(12) 



5 Noisy driving - Quantum Zeno effect 

Since quantum noise may be viewed as an average over all realizations of erratic 
driving, the expected decay in this case should follow Eq. p^ . We would like to 
obtain a more practical version of this formula using a slightly different procedure, 
that is inspired by the analysis of the Quantum Zeno effect [23 2 4. 25 26| i^. Previous 
analysis |28I29| lead to Fermi-Golden-rule type expression for the decay of coherence: 

2 

S = eM~rlAD^,t] , - [cot\26))]^ . (13) 

which is valid for to <C (l/wj), and restricted to short times such that Du,t <^ 1. 

Our analysis does not assume that Dj^t <ti 1, but we still assume that D^t <^ 
log(iV), such that S is not far from unity, and the folding of the wavepacket around 
the sphere can be ignored. During each squeezing interval tn the angular variance 
grows as {r'^)t+to = (r^)t [1 + 2 cot^ (269) sinh^(u» ,/<£>)], before being reset by the noise. 
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Fig. 3. The one-body coherence S is plotted as a function of time, for A'^ = 50 particles 
and u = 2. On the left: the dashed/yellow line is the drive-free decay. Single realizations of 
erratic driving whose intensity is D = 37.5wj are illustrated by gray lines. The red lines 
mark the 25%, the 50% (median), and the 75% percentiles over 2000 realizations. The dash- 
dotted/blue line is the average as obtained from the master equation Eq. ((S)). The symbols 
are based on Eq.[TT](o) and on Eq.[l2](A). On the right the theoretical predictions are tested 
for shorter times during which coherence is not yet lost. The parameters are A'' = 100 and 
u = 2. The noise-free decay is marked by the dash-dotted/black line, while the dashed/blue 
is for D = lOwj, and the solid/red is for D = AOwj. The symbols are based on Eq. [5] (A) 
and on Eq.[Tl(o, □). 

Approximating the sinh^ by a quadratic function, we get after t/tu steps an expo- 
nentially growing (r^)t, leading to 



The accuracy of this expression is demonstrated in the right panel of Fig. [31 



6 Summary 

We have studied the effect of driving on the hyperbolic instability of the Bose-Hubbard 
dimer. High frequency off-resonant harmonic driving results in the many-body equiv- 
alent of the Kapitza pendulum effect. On the other hand, erratic driving and quantum 
noise slow down the loss of single particle coherence via a many-body quantum Zcno 
effect. However, for long times the decay of S departs from the Fermi-golden-rule 
expectation. Namely, the interplay of angular diffusion with the hyperbolic squeez- 
ing results in log-wide statistics for S, whose typical value differs from the algebraic 
mean. 

This research was supported by the Israel Science Foundation (Grants No. 346/11 and No. 29/11) 
and by the United States-Israel Binational Science Foundation (BSF). 
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